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ABSTRACT 

We present initial results from a population synthesis model aimed at determining 
the star formation rate of the Milky- Way. We find that a total star formation rate of 0.68 
to 1.45 M Q /yr is able to reproduce the observed number of young stellar objects in the 
Spitzer /IRAC GLIMPSE survey of the Galactic plane, assuming simple prescriptions 
for the 3D Galactic distributions of YSOs and interstellar dust, and using model SEDs 
to predict the brightness and color of the synthetic YSOs at different wavelengths. This 
is the first Galaxy- wide measurement derived from pre-main-sequence objects them- 
selves, rather than global observables such as the total radio continuum, Ha, or FIR 
flux. The value obtained is slightly lower than, but generally consistent with previously 
determined values. We will extend this method in the future to fit the brightness, color, 
and angular distribution of YSOs, and simultaneously make use of multiple surveys, to 
place constraints on the input assumptions, and reduce uncertainties in the star for- 
mation rate estimate. Ultimately, this will be one of the most accurate methods for 
determining the Galactic star formation rate, as it makes use of stars of all masses (lim- 
ited only by sensitivity) rather than solely massive stars or indirect tracers of massive 
stars. 

Subject headings: infrared: stars — Galaxy: fundamental parameters — stars: forma- 
tion — stars: pre-main sequence 



1. Introduction 

Virtually all estimates of the current star formation rates (SFRs) for galaxies in the U niverse 



rely on stellar population synthesis models for calibration (as reviewed by iKennicuttl 1 19981 ) . These 



models predict global observables, such as the total Ha or infrared flux, that can be compared to 
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observations (e.g. ICalzetti et alj|2007l : iKennicutt et al.ll2009l ). This is necessary because individual 
forming stars cannot be resolved in any galaxy other than our own. 

Until recently, even estimates of the SFR of the Milky- Way have relied on global observables. 
Such studies generally rely on indirect tracers of massive (O- and early-B-type) stars to determine 
a massive SFR. This value is then extra polated to lower masses to derive a global SFR for the 
Galaxy. For example, ISmith et al.l (|1978l ) found a value of 5iWQyr , by making use of the fact 
that the integrated flux density from an H II region is a direct measure of the number of ionizing 
photons required to maintain that H II region, and is therefore an indirect measure of the number 
of O and e a rly-B -type stars. A more recent example includes the estimate of 4M yr _1 from 
Diehl et al.l ( 20061 ) who use a measurement of the total mass of 2.8 M & of 26 Al in the Galaxy 
from th e 7-ray flux mea s ured by INTEGRAL as a proxy for the massive star population of the 
Galaxy. iMisiriotis et al.l ( 20061 ) find a value of 2.7M0yr , by using the total 100 ^m flux of the 
Galaxy and a c onversion factor that depends o n the same population synt hesis model as used for 
other galaxies (|Misiriotis et al.l |2004j) . Finally, Murray &: Rahman! (J2009|) use the total free- free 
emission in the WMAP foreground map as a probe of the massive star population, and derive a 
global SFR of 1.3 M Q yr _1 . The difference between the values may in part be due to the fact that 
these studies do not all assume the same initial mass function (IMF). For a recent review of the 
tr acers and determinat ions of the Galactic SFR, normalized to the same IMF, we refer the reader 
to 



Calzetti et al.l (|2009T ). 



In general, any total flux measure that is related to the SFR of a galaxy (including the Milky- 
Way) is completely dominated by the high mass stars, since these are responsible for virtually all 
of the luminosity of a galaxy. Therefore, all of these methods are fundamentally limited by an 
extrapolation to lower masses, and higher sensitivity observations will not improve the estimates. 
For example, the ionizing flux for the H II regions or the 7-ray flux from 26 Al mentioned above are 
entirely dominated by the massive stars. 

In this letter we propose to measure the present-day Galactic SFR by directly observing and 
counting pre-main sequence stars. Estimates of the SFR of indi vidual regions with in a few hun- 
dred parsecs of the Sun have recently been derived in this way (jEvans et al.ll2009l ). yet no such 
measurement exists on the scale of the Galaxy. While Galactic plane surveys with finite sensitivity 
preferentially favor higher mass stars, as surveys reach higher and higher sensitivities, the IMF 
extrapolation to lower masses will gradually account for less and less of the resulting SFR, and the 
accuracy of this method will be much higher than for values derived from massive star formation 
indicators. In fact, this method may ultimately provide a more accurate measurement than for any 
other galaxy, because it is a direct measure of the very young, just formed stellar population, which 
is unresolved in most other galaxies. 
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Observations 



In this letter we m ake use of the Spitzer /IRAC GLIMPSE survey of the Galactic mid-plane 
(iBeniamin et alJ booah . specifically GLIMPSE I (covering 10° < \£\ < 65° and |6| < 1°) and 
GLIMPSE II (which fills in the region for \£\ < 10°, with \b\ < 1° for \£\ > 5°, \b\ < 1.5° for 2° 
< \£\ < 5°, and |6| < 2° for \£\ < 2°). The total area covere d by these two surveys is 274 deg 2 . We 
use the census of intrinsically red mid-infrared sources from iRobitaille et al.1 (120081 . hereafter R08) 
which includes 18,949 sources from these two surveys. 

The sources from R08 were selected using simple brightness and color selection criteria, namely 
13.89>[4.5]>6.50, 9.52>[8.0]> 4.01, and [4.5]-[8.0]>l. Extended sources are excluded. In total, 
18,949 sources were selected, and the completeness was estimated to be at least 65.7%. As shown 
in R08, these sources consist mostly of young stellar objects (YSOs; 50 — 70 %) and asymptotic 
giant branch (AGB) stars (30 — 50%), and may include a small number of planetary nebulae (at 
most 2 — 3%). Thus, there are likely to be approximately 9,500 to 13,300 YSOs in the R08 census. 
Taking into account the fact that the 65.7% completeness is a lower limit, the actual number of 
YSOs that would have been detected if the survey had been complete is in the range 9,500 to 
20,200. 



3. Model 

To derive a Galactic SFR using the R08 census, we construct a population synthesis model for 
YSOs in the Galaxy, and apply the same observational constraints as for the R08 census, varying 
the model SFR such that the number of 'detected' synthetic YSOs matches the observed number. 

The model is constructed as follows. First a 3D distribution for star formation is assumed, 
from which the YSO positions are randomly sampled. Each YSO is then assigned a random age 
and mass by assuming an IMF and a lower and upper stellar age. The number of synthetic YSOs 
sampled is thus directly related to the SFR for that model, which is given by the total mass of 
synthetic YSOs divided by the difference between the upper and lower age. Each synthetic YSO 
is then assigned intrinsic magnitudes at IRAC wavelengths (based on the stellar mass and age), 
which are scaled to its distance from the Earth. A 3D dust distribution and an extinction law are 
assumed in order to compute the extinction along the line of sight. Finally, the number of synthetic 
YSOs that fall inside the GLIMPSE I and II survey areas, satisfy the brightness and color selection 
criteria from R08, and are point sources, is compared to the observed number, and the SFR of the 
model is adjusted so that the two values match. The following paragraphs describe the assumptions 
made for each of these steps. 
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3.1. 3D YSO distribution 



The exact 3D distribution of YSOs in the Galaxy is unce rtain, but there exist stat istical 
estimates of where we expect stars to be forming. In particular, iBoissier Prantzosl (|1999l ) used 
the distribution of H and H2 as a function of Galactocentric radius to estimate the star formation 
rate relative to that at the solar radius, by assuming a Schmidt-type law relating the gas density 
and the star formation rate (SFRoc Xg a g), and taking into account the rate of passage of the gas 
through the spiral arms. The resulting radial profile was in good agreement with the distribution 
of tracers of massive stars (specifically H II regions, pulsars, and supernova remnants) found by 
previous studies. This radial function peaks around 4 to 5 kpc. T he values used were obtained from 
digitizing the solid line in Figure 2 of lBoissier &: Prantzosl (| 19991 ). 



We assume that the azimuthal distribution of YSOs is uniform, even though most star for- 
mation occurs in spiral arms. We make this assumption because the biggest effect of including 
spiral arms in the model would be to cause the longitude distribution of the sources to change on 
small scales, but would not change the overall number of detected synthetic sources significantly. 
Therefore, the azimuthally symmetric distribution we assume is still statistically representative of 
the overall distribution of YSOs. We assume a scale height of 80 pc, approximately the height of 
the Galactic thin disk. 



We take the Sun to be located a t a Galactocentric radius of 8.5 kpc (IKerr k, Lynden-Belllll986l ) 
and 27 pc above the Galactic plane (jChen et al.ll200ll ). 



3.2. Stellar masses and ages 



We adopt the IMF from iKroupal ( 200ll . equations (1) and (2)). Since the IMF is being used to 
implicitly extrapolate the observed mass to lower masses, it is also necessary to choose a minimum 
and maximum mass. The exact limits used are not important, as long as virtually all the observed 
YSOs fall inside that range. The SFR derived will then be the total SFR inside that range. We 
choose a lower mass of 0.1 Mq and an upper mass of 50 M , since we do not expect to see any 
objects outside these limits - 0.1 Mq sources would be too faint, and 50 M sources are very rare 
(and saturated in GLIMPSE). 

Recent observations suggest that primordial c ircumstellar disks are present for a few Myr 
then disappear on relatively short timescales (e.g. iFurlan et al.ll2009l ). Therefore, since we are 
considering YSOs with infrared excesses typical of primordial circumstellar disks (and of younger 
stages of evolution) we can assume that most of the R08 sources are younger than a few Myr. We 
choose a upper limit on the age of the synthetic YSOs of 2 Myr, and a lower limit on the age of 
1,000 yr, as we do not expect to be able to see any sources younger than this. It is unlikely that the 
overall star formation has changed significantly over the last few Myr, which is a small timescale 
compared to the Galactic rotation period of over 200 Myr. Therefore, we do not need to worry 
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about time variations in the SFR, in contrast to methods based on main-sequence stars, where 
choosing a constant SFR is a strong assumption. 



3.3. Dust 



The 3D distribution of dust in t he Galaxy is the subje ct of active research (e.g. ISale et al 



20091 ). We use the distribution from Misiriotis et al.1 (j200fil . Equations (1), (3), and Table (2)), 
which is a simple double exponential model based on modeling the far-IR dust emission, and is 
accu rate on large scales. There exist more complex models for the distribution of dust, such as that 
from lDrimmel Spergell (|200ll ) which includes spiral arms. The effect of a more clumpy extinction 
distribution will be to introduce scatter in the extinction values, but the overall effect is not likely 
to be significant. 

In recent years there have been a number of de terminations of the extinction law at mid- 
infrared wavelengths for both the interste llar medium (|Lutzlll999l : llndebetouw et al.ll2005l ) and for 
star forming regions ([Flaherty et al.ll2007l ). For this work, we are m odeling the large-scale ex tinction 
from the interstellar medium and therefore we use the values from llndebetouw et al.1 (|2005l ) for the 
IRAC bands. 



3.4. Intrinsic magnitudes of YSOs 

The final piece of the puzzle that is needed to construct our model is a prescription to convert 
the masses and ages of the sources into observable magnitudes. Unfortunately, no simple prescrip- 
tion exists for this, as the brightness and color depend on the amount of circumstellar dust and the 
properties of the central source, which in turn are far from certain for all stellar masses and ages. 
Of all the assumptions made in this letter, this is likely to be most crucial one. 



In lRobitaille et al.l (|2006l . hereafter R06) we computed a large set of model SEDs for YSOs with 
a large range of evolutionary stages (from heavily embedded sources to transitional disks) and stellar 
masses (from 0.1 M to 50 M ). To sample the parameters, each source was assigned a mass and 
age, and the stellar temperature and radius were found from evolutionary tracks. Subsequently, 
the disk and envelope parameters were sampled as a function of these stellar parameters. The 
sampling attempted to take into account observational and theoretical constraints on pre-main- 
sequence evolution. It is important to note that the parameter space in R06 does not constitute an 
evolutionary scenario, but rather allows for different scenarios, some of which may be unrealistic. 
However, the general trends are nevertheless close to how we believe stars form. 

For each synthetic YSO, we choose the model with the closest stellar mass and age, and sample 
a random viewing angle. Since we are assuming that the sources in the R08 census are younger than 
2 Myr, the R06 models with ages between 2 and 10 Myr are not used. The IRAC magnitudes of the 
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models we do use are dominated by the stellar luminosity and temperature, and the envelope infall 
rate (which determines the envelope mass and therefore the degree of obscuration of the central 
star). Therefore, the most important assumptions from R06 that will impact this work are the 
choice of stellar evolutionary tracks (which fix the stellar luminosity and temperature for a given 
age and mass) and the choice of the time dependence of the envelope infall rate. The latter is 
sampled from a range that is constant before 0.1 Myr and decreases to zero by 1 Myr. 

The R06 models were computed in 50 circular apertures from 100 to 100,000 AU, effectively 
providing radial brightness profiles at all wavelengths. Therefore, in the population synthesis model, 
we can find the flux for each synthetic source in an aperture corresponding to the resolution of the 
IRAC observations (2"). In addition, we use this information to reject sources that would have 
been extended in the GLIMPSE survey, and would therefore have been excluded from the R08 
census. We do this by selecting only sources where 99 % of the flux falls inside an aperture 2" in 
radius. This reduced the number of final 'detected' YSOs by approximately 10 %. 

4. Results 

Ignoring uncertainties in the input assumption, we find that a SFR in the range 0.68 to 
1.45 M /yr reproduces the number of observed YSOs, where the range of values accounts for the 
uncertainty in the actual number of YSOs, which is in the range 9,500 to 20,200 (^}2]). To compare 
the properties of the 'detected' synthetic YSOs with those of the 11,649 selected candidate YSOs 
from R08, we computed a model with 2.73 x 10 6 synthetic YSOs - corresponding to a SFR of 
0.83 M Q /yr - of which 11,919 were 'detected'. This model is the one shown in Figures Q] to [3j 
Figure Q] shows the mass distribution of all the synthetic YSOs, the subset that fall in the survey 
area, and the subset that would have been detected by Spitzer and included in the R08 census. 
The distribution of 'detected' synthetic sources is biased towards high masses, and peaks just under 
10 Mq. Around 50% of all 10 M Q synthetic YSOs are detected. The drop-off at lower and higher 
masses is due to the lower and upper brightness cutoffs. Since approximately 2.7 million synthetic 
YSOs between 0.1 and 50 M are required Galaxy- wide to explain the observed numbers, the R08 
census represents less than 0.5% of all YSOs in the Galaxy. 

Figure [2] shows the 3D spatial distribution of all the synthetic YSOs that would have been 
included in the R08 census, color-coded by the interstellar extinction to Earth. Most sources 
are seen within lOkpc but many are seen out to the far side of the Galaxy, even though the 
interstellar extinction is significant. The ring with a 5kpc radius is due to the peak in the radial 
SFR distribution assumed ( §3.1[) . but we note that the real distribution of YSOs does not necessarily 
trace a ring. 

Finally, Figure [3] shows the longitude and latitude distribution, as well as the [8.0] and [4.5] — 
[8.0] distributions for all the synthetic YSOs that would have been included in the R08 census, and 
for all the YSO candidates in the R08 census. Even though no effort was made to fit the model 
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distributions to the observed ones, the two are in reasonable agreement. The [8.0] magnitude 
distributions do diverge at the faint end, but the difference never exceeds 50 %, and this difference 
may be due to remaining contaminating AGB stars. Both the observed and model longitude 
distribution drop off around longitudes of 30 to 40° on either side of the Galactic center. These 
longitudes correspond to those where the line of sight is tangent to the peak in the radial distribution 
of star formation ( §3.ip . Also of interest is that the observed latitude distribution of sources is 
asymmetric, with 55.9 ± 0.6 of sources at negative latitudes. The model also predicts a similar 
asymmetry (albeit slightly larger), with 61.5 ± 0.7% 1 of sources at negative latitudes. The reason 
that the model presents an asymmetry is entirely due to the fact that the Sun is displaced 27 pc 
above the Galactic mid-plane. Therefore, it is safe to infer from this that the asymmetry in the 
latitude distribution of observed sources is also due to the displacement of the Sun above the 
Galactic mid-plane. 



5. Summary 



By developing a population synthesis model for the Milky-Way, and adjusting the star for- 
mation rate such that the number of synthetic YSOs that would have been present in the R08 
census matched the actual number of observed YSOs, we obtain a star formation rate of 0.68 to 
1.45M©/yr. The uncertainties on this value take into account only the uncertainty in the actual 
number of observed YSOs (since these are not perfectly separable from AGB stars) and the uncer- 
tainty in the completeness. Our result is of the same order as but lower than previous estimates 
(©> but this can be explai ned in pa r t by diff erences in th e assumed IMF. The previous values 



quoted in 32 assume either a lSalpeterl (|1955l ) or IScalol (|1986l ) IMF, whi ch over-predict the number 
of low-mass stars by approximately 50% compared to a more realistic iKroupal ( 200ll ) IMF. If we 



assume a Salpeter IMF from 0.1 to 50 M©, we obtain a SFR in the range 0.98 to 2.09M©/yr, in 
better agreement with previous values. We find that the YSOs in the R08 census are dominated by 
sources in the range 3 to 20 M®, and that these are seen out to very large distances (10 to 15 kpc) in 
some cases. The longitude and latitude distribution of the synthetic YSOs are in reasonable agree- 
ment with the observed distribution, and we find that the asymmetry in the latitude distribution 
can be entirely explained by the fact that the Sun is not located exactly in the mid-plane. 

The input assumptions for the model are simple, consisting of an axisymmetric distribution of 
synthetic YSOs peaking around 4 to 5 kpc, a double exponential distribution for the dust, a Kroupa 
IMF, a standard mid-infrared extinction law, and a prescription for the magnitudes of YSOs based 
on the stellar mass and age, using the models from R06. The assumption that is likely to matter 
the most is the latter, namely knowing what brightness and color to assign to YSOs as a function of 
age and mass. This is because while there are many examples of well-studied low-mass YSOs with 
approximate ages, the pre-main-sequence evolution of 10 M© stars, which dominate the 'detected' 



1 where the uncertainty is from Poisson statistics 
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synthetic YSOs, is poorly determined. 

The model will be developed into a much more powerful diagnostic tool in the future. In 
particular, we plan to fit the angular, brightness, and color distributions of the sources to improve 
the estimate of the SFR and place constraints on other input parameters and assumptions. For 
example, the angular distribution of the sources will help place constraints on the model for the 
SFR as a function of Galactocentric radius, and the color and brightness distributions will help 
understand whether the assumptions for the intrinsic magnitudes of YSOs and the extinction 
properties are sensible, or whether they need to be modified. In addition, we will thoroughly 
explore the parameter space in order to understand how well constrained the input assumptions 
and parameters are, and what additional data could resolve degeneracies. Finally, we plan to model 
multiple datasets simultaneously, to reproduce not only the observable properties of sources seen 
in any one survey, but all available infrared or sub-mm surveys of the Galactic plane, such as the 
UKIDSS and Herschel surveys. 
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Fig. 1. — The mass function of the synthetic YSOs for one particular run of the population synthesis 
model. The yellow histogram shows all of the synthetic YSOs, the orange histogram shows only 
the sources that fall in the survey area, and the red shows only the sources that would have been 
detected by Spitzer and included in the R08 census 
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Fig. 2. — The 3-D spatial distribution of the synthetic YSOs that would have been detected by 
Spitzer and included in the R08 census, for one particular run of the population synthesis model, 
color-coded by the interstellar Ay to the Earth (c.f. colorbar) 
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Fig. 3. — The Galactic longitude (top left), Galactic latitude (bottom left), [8.0] magnitude (top 
right), and [4.5] — [8.0] color (bottom right) distribution for the observed YSOs from the R08 
census (gray filled histogram) and for the synthetic YSOs from one particular run of the population 
synthesis model (red solid line histogram). The dashed line in the Galactic latitude panel shows 
b = 0° to emphasize the latitude asymmetry. The dashed lines in the magnitude and color panels 
show the selection criteria used in R08. 



